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We analyze an interplay between Coulomb blockade and quantum fluctuations in a coherent 
conductor (with dimensionless conductance g ^ 1) attached to an Ohmic shunt. We demonstrate 
that at r = the system can be either an insulator or a metal depending on whether its total 
resistance is larger or smaller than h/e'^ ~ 25.8 kft. In a metallic phase the Coulomb gap is fully 
suppressed by quantum fluctuations. We briefly discuss possible relation of this effect to recent 
experiments indicating the presence of a metal-insulator phase transition in 2d disordered systems. 



It is well known that Coulomb interaction may strongly 
affect quantum transport of electrons in mesoscopic con- 
ductors. For instance, electron tunneling across metal- 
lic junctions can be strongly suppressed or even blocked 
completely at T = due to Coulomb effects [|l| ^. This 
Coulomb blockade of tunneling is a direct consequence of 
the electron charge discreteness, manifestations of which 
persist even for junctions with resistances well below the 
quantum resistance unit Rn — h/e^ ~ 25.8 kl7. 

Recently it was argued [Q that the I — V curve of an 
arbitrary - albeit relatively short - coherent conductor 
in the presence of interactions can be expressed in the 
form 



R^ = l-PfiV,T), 



(1) 



where 1/R — (2e^//i) T„ is the Landauer conduc- 
tance of a scatterer and T„ are the transmissions of its 
conducting modes. The magnitude of the interaction 
term in eq. (bj) is controlled by the parameter 



E„rn(i-r») 
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and f{V,T) is a universal function which depends on R 
as well as on the external impedance Zs{uj) (e.g. leads) 
attached to the scatterer. A similar result was also ob- 
tained in Ref. |^ in the limit of a single conducting mode, 
in which case (3 — 1 — Ti . 

The parameter (|^) is already well known in the theory 
of shot noise Q . Since shot noise and interaction effects 
in mesoscopic conductors arc both the manifestations of 
the electron charge discreteness, it appears quite natural 
that the magnitude of both these effects is governed by 
the same parameter j3 (|^). 

It is also important that the above effects do not dis- 
appear even if the charge quantization is not preserved 
in a part of the system (e.g. in the leads). Consider 
an important example of an Ohmic external impedance 
Zs{u>) ~ Rs and define gs = Rq/Rs and go = g + gs. 
The function f{V,T) was already evaluated before for 



the case of tunnel junctions. At not very high energies 
max(T, eV) ^ goEc and for go ^ 1 oire finds 



f{V,T) 2, {2/go)\n{goEc/max{T,eV)), 



(3) 



i.e. in this regime Coulomb interaction causes a loga- 
rithmic correction to the scatterer conductance. This 
correction is small in the parameter l/go- Here Ec is the 
effective charging energy which sets the energy scale for 
the interaction effects in our problem. The result (^ is 
valid down to energies max(r, eT^) ^ gg^c 6xp(— (7o/2). 
At even lower T and eV energy relaxation effects turn 
out to be particularly important making the conductance 
saturate at a universal value B 



G = 



I- 13 _ 2e2 
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This formula recovers complete Coulomb blockade in the 
limit (3 ^ \ (tunnel junctions), demonstrates the absence 
of it for ballistic scatterers (/? — > 0), and yields suppres- 
sion of the Landauer conductance by the factor 2/3 in 
the important case of diffusive conductors. 

Within the approach J |] one can demonstrate that 
all corrections to eq. (^) are small in the parameter 
(3 /go <C 1. Nevertheless, one cannot yet conclude that 
the system behavior is metallic at T = because the 
Langevin equation analysis does not fully account for 
the charge discreteness in our system. Specifically, this 
analysis does not include subtle nonperturbative effects 
which can be captured only by means of the instanton 
technique. These effects may provide a mechanism which 
turns a conductor into an insulator at T = 0. Thus, in 
order to understand the ground state properties of a sys- 
tem with disorder and interactions, the instanton effects 
should be properly taken into account. 

In this Letter we will analyze the role of the instan- 
ton effects in relatively short coherent conductors. An 
important part of the work was already carried out by 
Nazarov Q who found the renormalized Coulomb gap 
for a coherent conductor in the form Ec oc Ec exp(— ag), 
where a = 7r^/8 for diffusive conductors and a = 1/2 
for tunnel junctions H. At temperatures below Ec the 
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charge quantization plays a crucial role and should yield 
insulating behavior at T = 0. We will extend the in- 
stanton analysis |Q in several respects. Our main goal 
is to study the effect of an external impedance Zg on 
the renormalized Coulomb gap Ec- In particular, we 
will demonstrate that for < Rq quantum fluctua- 
tions yield complete suppression of the Coulomb gap Ec 
and, hence, destroy the insulating state even at T = 0. 

The model and effective action. As in Ref [Q we will 
consider an arbitrary coherent scatterer between two big 
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reservoirs. Phase and energy relaxation are only allowed 
in the reservoirs and not during scattering, i.e. the scat- 
terer is assumed to be shorter than both dephasing and 
inelastic relaxation lengths. Coulomb effects in the scat- 
terer region are described by an effective capacitance 
C. The charging energy Ec — e^/2C is assumed to be 
smaller than the typical inverse traversal time. 

The grand partition function Z of the system "scat- 
terer -|-environment" and the scatterer effective action 
can be written in the form |^,|0 11 
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sin[7rr(T2m - Ti] 
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Here <p(T)/e = V[t) and q = I{t) are respectively 
the fluctuating voltage across the scatterer and the cur- 
rent through the system. Eq. was initially derived 
for perfectly transparent contacts JTof and generalized in 
Ref. pTl to arbitrary transparencies. One has 
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where A is the cross section of the contact and D{0) is 
its angle-dependent transparency. The sum ^.^^ is taken 
over independent conducting channels and i?„ ~ 1 — r„. 

Isolated scatterer. Let us first consider an isolated scat- 
terer 1 /Zs and rederive the result [Q for Ec from our 
effective action (^. In this limit integration over q in (|^) 
fixes the charge to be constant i J dT{ipq/e) 2i:ilqx/e. 
At large conductances 5^1 the remaining integral over 
Lp is evaluated within the saddle point approximation. 
We are interested in the nontrivial saddle points (p{t) for 
the action which "connect" states with different wind- 
ing numbers A general expression for such instantons 
is given in Q and in a certain limit it reduces to a set of 
Korshunov's instantons ||l^ i^(r) = 2arctan(n(T — tq)). 
Since the saddle point action S'o[i^] does not depend on il, 
it suffices to set ^ 0. Then (p reduces further to a set 
of "straight hnes" 2'kTtI Substituting (^(r) = 2ttTt 
into (0) one trivially gets 
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Combining this equation with Ec /Ec oc exp(— S'o[(^(t)]) 
one arrives at the result (for a spin degenerate case). 

It is also possible to go beyond the exponential accu- 
racy and to estimate the pre-exponent in the expression 
for Ec. Here we restrict ourselves to an approximation 



of non-interacting instantons. Within this approximation 
Ec can be found by means of a simple formula |14 



N 



T 



(5o[^(r)]))^/2exp(-5o[^(r)]), (9) 



which yields correct results up to an unimportant numer- 
ical prefactor of order one. Here Lk and Zk are respec- 
tively the effective volume and the effective instanton size 
for the fc-th zero mode and N is the total number of zero 
modes. Similarly to in our problem each instanton 
has two zero modes, corresponding to shifts of its center 
To in time from zero to 1 /T and to fluctuations of its fre- 
quency r2 within the interval < 17 ^ Ec (fluctuations 
with Vl > Ec are exponentially suppressed by the charg- 
ing term in the action). Thus we have N = 2, Li = 1/T 
and L2 ~ Ec. The parameters zi^2 for both zero modes 
can be easily evaluated. Substituting the eigenfunctions 
for these zero modes dcp/dr and dip /d ft into the effective 
action one finds zi ~ 1/Q and Z2 ~ Thus, the prod- 
uct Z1Z2 is ri-independent and just reduces to a numerical 
factor of order one Z1Z2 1. Then eq. (^) immediately 
yields 



Ec/Ec 



In 



ag exp(~ag). (10) 



This formula is valid for ag 3> 1, i.e. either at large con- 
ductances g 3> 1 or, if g ^ 1, for very small values i?„ 
implying a ^ 1. In a (spin-degenerate) single channel 
limit eq. (|o|) reduces to Ec/Ec RilnR^^ in agree- 
ment with the result derived in Ref. |l^] for i?i <C 1 
within a different technique. 

In order to find the dependence of the ground state 
energy Eq on q^ it is necessary to sum over all possible 
instanton configurations. For non-interacting instantons 
this summation yields 
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In the particular case of a tunnel junction in the strong 
tunneling limit 17 ^ 1 this result agrees with one de- 
rived in eq. (10) of Ref. §, where the pre-exponent in 
the expression for A (simply proportional to the inverse 
i?C-time) was obtained by means of an explicit calcula- 
tion of the fluctuation determinants. We also note that 
(comparatively weak) inter-instanton interaction might 
slightly modify the result (pi]), both the form of Eo{qx) 
and the pre-exponent in the expression for A. However 
this effect will not be important for us here. 

Now we are in a position to study the scatterer con- 
ductance at extremely low energies. In the limit of high 
external impedances Zs 3> Rq both the resistance i?(w) 
and conductance G{uj) — l/i?(w) of our scatterer are de- 
termined from the correlation function {fip) evaluated 
for 1/Zs 0. The calculation is again performed in 
imaginary time and is completely straightforward. Simi- 
larly to Ref. [0 let us make a shift Qx ^ qx+ £,{t) in the 
exponent of eq. (|^) and define the phase-phase correlator 
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(12) 



In order to evaluate the generating functional Z[^(^)] (|^) 
for ag ^ 1 it is sufficient to consider small fluctuations 
of if and instantons. Since we are mainly interested in 
the low frequency behavior of the correlator (^) we can 
safely assume that typical frequencies of the source field 
do not exceed the instanton ones fl. Under this adi- 
abaticity condition our previous instanton analysis can 
be trivially repeated. Then for — > from eq. (12) we 
obtain 



e^Reff 47r2A 
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Here uj is the Matsubara frequency. What remains 
is to define the "Matsubara resistance" i?(w) — 
[\uj\/ e'^){(p{uj)ip{—uj)) and perform an analytic continu- 
ation to real frequencies. After that at T = one gets 



G(u;) 
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Reff + (27r/e)2A/iw' 
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This result implies that - while at 3> A the con- 
ductance of our system remains finite - in the opposite 
low frequency limit w ^ A the renormalized Coulomb 
gap A ~ Ec becomes important, the response is pre- 
dominantly capacitive and G(ijj) vanishes at — > 0. 
At nonzero but low T <^C A and 1/Zs one has 
G{T) oc exp(— A/r), i.e. the scatterer response is in- 
sulating at T = due to Coulomb blockade. A simple 
estimate of the effective resistance Reff in (p^, (14), 
similarly to |^ , would give Reff — R- In the spirit of 
our Langevin equation analysis B, one can also expect 
a (frequency dependent) correction to this estimate. In 



particular, at low u < gEc Gxj>{—g/2) it is natural to 
expect Reff — R/{^ — P)- 

Metal- Insulator phase transition. Now let us see how 
the above behavior is modified in the presence of a fi- 
nite external impedance Zs{oj). In this case the charge 
q fluctuates and, hence, should be treated as a quantum 
variable. For our purposes it is convenient to first per- 
form the path integral (|^) over the phase ip{t). As before, 
in the limit ag ^ \ we integrate over small fluctuations 
of Lp and instantons. Integrating out small fluctuations, 
for \uj\ ^ Ec one arrives at an effective impedance 



Z{-iu)^Zs{~'iuj) + R, 



(15) 



i.e. the scatterer impedance should simply be added to 
one of the shunt. Evaluating the instanton contribu- 
tion we again assume that typical frequencies uiq for the 
charge variable q do not exceed J7. This is sufficient for 
Zs 3> R- Then, similarly to P], we find 



Z = j P(jexp|^-|^|a;|Z(-zc^)|g^|2- j 



drEoiq) 



Let us now choose Zg to be Ohmic Zg{u!) ~ Rs- In this 
case the resulting effective action for the charge q coin- 
cides with one for a linearly damped quantum particle 
in a periodic potential (pT|). This is a well-known prob- 
lem [0,0 which can be treated, e.g., by means of the 
standard renormalization group (RG) technique. Suc- 
cessively reducing the high frequency cutoff tOc and inte- 
grating out charges with higher ujq one arrives at the RG 
equations Iit) for A = A/ujc and — 993 /go'- 



dA 



d(ln LUc) 



= A(.9s - 1) + 0(A3 



dgi 



d{ln ujc 



= 0. (16) 



Eqs. ( p^ ) demonstrate that for gs < 1 the value A - 
though initially small - grows in the course of renormal- 
ization. Thus the Coulomb gap remains finite and the 
system is an insulator at T = 0. On the other hand, 
for gs > 1 the Coulomb gap scales to zero. In this case 
Coulomb blockade is completely destroyed by quantum 
fluctuations and the system behavior should be metallic 
down to T = 0. 

The charge-charge correlation function can also be de- 
rived in a straightforward manner. At |ci;| < Ec one 
finds 



{qq)^ ~ ((i? + Rs)\oj\ + vr^A./e^)-!. 



(17) 



Here A^ plays the role of the renormalized Coulomb gap. 

It is equal to A^ = A(A/(jJc) ^""^ for .gs < 1 and A^ = 
otherwise. As expected, charges are localized in the in- 
sulating phase and delocalized in the metallic one. 

Choosing the initial cutoff frequency as uJco Ec, one 
can rewrite the first eq. ( p^ ) directly for the combination 
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d{ag)/d{lnLU^) = (1 - gj:){l + 1/ag). 



(18) 



This equation is valid as long as ag ^ 1. In the tunnel- 
ing limit T„ < 1 one has EJ^K^ - En^n = 9/2, 
and ([l8| ) reduces to the RG equation for the tunneling 
conductance derived in |^ as well as in Refs. [p"8|jl9| ] (for 
the case = 0). 

Combining our present results with those of Ref. Q| we 
arrive at the following picture. In the limit of large con- 
ductances go ^ 1 the I—V curve of an arbitrary coherent 
scatterer is described by eqs. (0)-(^ down to exponen- 
tially low temperatures and voltages max(T, ey) ^ A^. 
In this regime Coulomb effects lead to partial suppres- 
sion of the scatterer conductance which saturates at the 
value (^ at low T and V. For gs > 1 or if at least one 
of the channels is ballistic, _R„ = 0, the Coulomb gap is 
fully suppressed = 0, and the conductance remains 
nonzero (Q) even at T = 0. This is the metallic phase. 
On the other hand, if < 1 and all i?„ > 0, at energies 



max{T, eV) < ^ Ec[ag e:cp{—ag)] 
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Coulomb interaction leads to further suppression of the 
scatterer conductance G which eventually vanishes at 
T = 0. This is the insulating phase. A quantum 
phase transition between these two phases occurs either 
at gs = 1 or if the conductance (per spin) of at least one 
of the channels gn = Tn becomes exactly equal to one. 

Summarizing, an interplay between charge discrete- 
ness, coherent scattering and Coulomb interaction yields 
two types of effects. One effect is controlled by the pa- 
rameter (3 (|^) and results in partial suppression of the 
system conductance at low T. For go ^ 1 this sup- 
pression is fully captured by the quasiclassical Langevin 
equation approach Q . This effect of the electron-electron 
interaction is relatively robust in the sense that the in- 
teraction correction (3f (|l|) is never suppressed unless all 
the channels are ballistic (/? = 0) and/or the shunt con- 
ductance gs is infinite. Another effect is due to instan- 
tons which give rise to the renormalized Coulomb gap 
Ar- This interaction effect is much more subtle since the 
Coulomb gap can be destroyed by quantum fluctuations 
much easier. It is this suppression which yields a quan- 
tum insulator-to-metal phase transition in our system. 

Finally, let us briefly address possible implications of 
our results for recent experiments [^0| which strongly 
indicate the presence of a metal-insulator phase tran- 
sition in various 2d disordered systems. One can con- 
sider a (sufficiently small) coherent scatterer with the 
dimensionless conductance g viewing all other scatter- 
ers in the system as an effective environment with the 
conductance gs- On a phenomenological level one can 
assume this environment to be Ohmic at sufficiently low 
frequencies. Under this assumption one immediately ar- 
rives at the conclusion about the presence of a quantum 
metal-insulator phase transition at 5s=l. In 2d systems 
one has g ^ gs ^ 90 ^ Qt.- Therefore in such systems 



this phase transition should be expected at conductances 
~ 1/Rq, exactly as it was observed in numerous ex- 
periments Local properties of the insulating and 
metallic phases are expected to be very different. In the 
insulating phase charges should be localized around in- 
homogeneities (puddles) due to Coulomb blockade, while 
in the metallic phase the Coulomb gap is destroyed by 
quantum fluctuations and the charge distribution should 
be much more uniform. These expectations are fully con- 
sistent with recent experimental observations [ pl| . Thus, 
there might be a direct relation between the experimental 
results (2^j2^ and the old problem Q describing dissipa- 
tive dynamics of a quantum particle in a periodic poten- 
tial. In order to explore this possibility a more detailed 
theory, which would include, e.g., the issue of quantum 
decoherence, would be warranted. 
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